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ABSTRACT 


Porosity of functionally graded materials (FGMs) is usually aroused by fabrication defects. It had been proven 
that the porosity has a significant influence on the static responses of their structures, but the effects of porosity 
on buckling behaviors are still worth investigating. To reveal these effects, the thermal-mechanical coupling 
buckling issue of a clamped-clamped porous FGM sandwich beam is investigated in this paper by employing the 
high-order sinusoidal shear deformation theory. The modified Voigt mixture rule is used to approximate the 
temperature-dependent material properties of porous FGMs. The physical neutral plane of FGM sandwich beams 
is taken into account to reflect the actual condition of the structures and simplify the calculation. The thermal 
environments are considered as uniform, linear and nonlinear temperature rises, and both the temperature- 
independent and temperature-dependent material properties are discussed in order to justify the importance of 
the thermal-mechanical coupling effect. An iterative algorithm is used to solve the thermal-mechanical coupling 
critical buckling temperature. The present theoretical results are verified by comparing with the literature and 
ABAQUS results, and the effects of porosity, the physical neutral plane, gradient index, material temperature 
dependence, sandwich structural parameters are discussed. Results show that for buckling issue excluding the 
pre-buckling deformation effect, considering either the physical neutral plane or the geometrical middle plane of 
FGM beams would produce alike critical buckling temperatures. With the rise of porosity, the critical tem- 
perature increases greatly, which is quite different from the changing rule observed in the buckling issue of 
inplane-loaded porous FGM plates in literature. The beam with a smaller face-to-core ratio is more sensitive to 
the change in porosity. Moreover, to improve the thermal buckling load of FGM beams, ceramic constituents 
with the lower thermal expansion coefficient would be preferred. 


1. Introduction 


et al. [2] studied the thermal buckling response of functionally graded 
plates subjected to uniform, linear and nonlinear temperature rises 


Functionally graded materials (FGMs) are a new kind of composites 
fabricated by mixing ceramic and metal constituents. As ideal heat-re- 
sistant materials, FGMs possess both the thermal resistance of the 
ceramic constituent and the crack-resistance of the metal constituent 
simultaneously. Controlling the volume fraction of the constituents in 
space domain would achieve their smooth gradient material properties, 
and thus excite their potential applications in optimizing thermo- 
dynamic performances of their structures, such as high-temperature 
and high-stress components in aircraft and automotive engineering 
fields [1]. 

Many studies focused on the thermal-mechanical coupling stability 
issues of FGM structures in the recent years. For example, Bouhadra 
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across the thickness direction. They found that for all the three thermal 
loading types, the thermal force resultant becomes larger as the ratio of 
metal-to-ceramic modules increases. Trinh et al. [3] investigated the 
thermal buckling of FGM beams under various thermal loads based on 
the state space approach. The results revealed that the critical tem- 
peratures under nonlinear temperature rise case are larger than those 
under linear temperature rise case, and the difference becomes more 
pronounced as the volume fraction index of ceramics increases. Besides, 
the solution considering the temperature-dependence of material 
properties gives significantly lower critical buckling temperatures than 
those with temperature-independence. Using a two-step perturbation 
method, She et al. [4] studied the thermal buckling and postbuckling of 
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clamped-clamped FGM beams subjected to uniform temperature rise. 
The results indicated that the thermal post-buckling equilibrium paths 
for FGM beams are stable, and the thermal post-buckling response of 
FGM beams under temperature-dependent case is lower than that under 
temperature-independent case. 

However, due to the great difference in the solidification tempera- 
tures between the ceramic and metal constituents, FGMs may have 
some micro-voids during the sintering preparation [5]. Wattanasa- 
kulpong et al. [6] also found the similar porous phenomenon of FGMs in 
multi-step continuous infiltration process. Investigations on the effects 
of porosity were addressed mainly for vibration characteristics of 
porous FGM beams [7-11], in which the modified Voigt mixing rule 
was used to approximate the material properties. The results indicated 
that porosity may lead to a great decrease in vibration frequencies. For 
buckling issues, Atmane et al. [12] and Barati et al. [13] presented the 
buckling analysis for in-plane compressed porous FGM beams and 
plates resting on an elastic foundation. It was proven that the critical 
buckling load is inversely proportional to the value of porosity, and 
considering the uniform porosity distribution in FGMs would give a 
smaller critical buckling load than the non-uniform distribution. 

Besides, most of the aforementioned literature established their 
theories based on the geometrical middle plane. But due to the asym- 
metry of material properties of FGMs with respect to the geometrical 
middle plane, the stretching and bending equations are coupled. With 
the physical neutral plane concept, Zhang [14] and Abrate [15] have 
shown that if the reference plane is properly selected, there is no 
stretching-bending coupling in constitutive equations. Prakash et al. 
[16] studied the effect considering physical neutral plane on the non- 
linear stability behavior of FGM plates under compression. It is found 
that the post-buckling path of FGM plates would be quite different 
when the physical neutral plane is taken into account. However, with 
an efficient sinusoidal shear deformation theory and the physical neu- 
tral plane concept, Bouiadjra et al. [17] investigated the thermal 
buckling behavior of FGM plates and the results have a good agreement 
with those based on the geometrical middle plane. Also, Bellifa et al. 
[18] presented bending and free vibration analysis of FGM plates using 
a simple first-order shear deformation theory based on the exact posi- 
tion of the physical neutral plane, and the results are identical to those 
of the conventional first-order shear deformation theory. Moreover, 
more and more researches have been attached to investigate the effects 
of the physical neutral plane in bending, free vibration, nonlinear vi- 
bration issues of FGM structures in the recent years [19-25]. 

To the author's best knowledge, there is no investigation on the 
thermal-mechanical coupling buckling behaviors of porous FGM beams, 
especially for porous FGM sandwich beams. In this paper, thermal- 
mechanical coupling buckling issue of clamped-clamped porous FGM 
sandwich beams is investigated by employing the high-order sinusoidal 
shear deformation theory (SSDT). The modified Voigt mixture rule with 
porosity is used to approximate the temperature-dependence material 
properties of porous FGMs and the physical neutral plane of FGM beams 
is taken into account. The thermal environments are considered being 
uniform, linear and nonlinear temperature rises. Both temperature-in- 
dependent and temperature-dependent material properties are ex- 
amined. Finally, an iterative algorithm is used to solve the thermal- 
mechanical coupling critical buckling temperature. 


2. Porous FGM sandwich beams 


Consider a clamped-clamped porous FGM beam with length L and 
width b and height h = he + 2h, as shown in Fig. 1. It is supposed that 
the displacement of the beam between the layers is continuous. The top 
and bottom layer of the sandwich beam are respectively made of 
homogeneous ceramic and metal constituents with the height h,, while 
the core layer is made of FGMs with the height hy. By setting k = 1,2,3 to 
denote the top, core and bottom layer of FGM beams, respectively. And 
we have the volume fraction of ceramic constituents V“ in the FGM 
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Fig. 1. Geometry of a clamped-clamped porous FGM sandwich beam. 
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where the subscript cm denote the ceramic and metal constituents 
respectively and n is the power-law index indicating the inhomogeneity 
of FGMs. 

Additionally, according to the principle of the multi-step sequential 
infiltration technique that can be employed to fabricate FGM samples 
[6], the porosities mostly occur at the middle zone due to different 
infiltration process. Thus, two types of porosity distributions in core 
named as FGM-I and FGM-II are discussed. For FGM-I, the porosity 
phases distribute uniformly over the cross section, while for FGM-II, the 
porosity phase distributes around the mid-zone of the cross-section in a 
rhombus shape, as shown in Fig. 1 (b). With the modified Voigt mixing 
rule proposed by Wattanasakulpong et al. [8], the effective material 
properties of FGMs P(z) representing the Young's modulus E(z), 
thermal expansion coefficient a(z) and thermal conductivity x(z) can 
be approximated as 


PQ (z) = (P. — OG + 2) + Pa- ce + Pa) A 


where £ (6 « 1) denotes the volume fraction of porosity, distributed 
equally in the ceramic and metal phases. And we have ¢ = 1 for FGM-I 
and ¢ = 1 — 2lzl/h~ for FGM-II. 

It is clear that the effective material properties decrease with the 
increase of porosity for both FGM-I and FGM-II. Meanwhile, the ma- 
terial properties should be temperature-dependent and can be ex- 
pressed as 


P, = co(c1Tı + 14+ G&T + (oT? + T?) (3) 


where i = c, m. T is temperature in Kelvin, and c;(j = —1,0,1,2,3) are the 
material temperature coefficients [3,26] listed in Table 1. The tem- 
perature-dependent material properties are also illustrated in Fig. 2. 
The Poisson's ratio of FGMs is assumed to be constant, i.e. 
V = (v% + Ym)/2 evaluated at T = 300K. 


3. Theory and formulation 
3.1. Basic formulations 


The position of the physical neutral plane (z = Co) satisfies the 
following relation. 
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Table 1 
Temperature-dependent coefficients for ceramic and metals [3,26]. 
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Constituents 


Properties Co C1 c2 C3 
SisN4 E(Pa) 348.43 x 10° — 3.070 x 1074 2.160 x 10-7 — 8.946 x 107! 
a(1/K) 5.8723 x 1076 9.095 x 1074 0 0 
x(W/mK) 13.723 = 1.032 x 1073 5.466 x 1077 — 7.876 x 107" 
v 0.2400 0 0 0 
ZrO2 E(Pa) 244.27 x 10° — 1371 x 10-3 1.214 x 10-6 — 3.681 x 10710 
a(1/K) 12.766 x 1076 — 1.491 x 10-3 1.006 x 1075 — 6.778 x 107! 
x (W/mK) 1.7 1.276 x 1074 6.648 x 1078 0 
ù 0.2882 1.133 x 1074 0 0 
SUS304 E (Pa) 201.04 x 10° 3.079 x 1074 — 6.534 x 1077 0 
a(1/K) 12.330 x 1076 8.086 x 1074 0 0 
x(W/mK) 15.379 — 1.264 x 1073 2.092 x 1076 — 7.223 x 10710 
v 0.3262 — 2.002 x 1074 3.797 x 10-7 0 
5 , 
Ux (x, Z) = u(x) — Z — Co)wo(x) + [P(Z) - C] po~), uz, z) 
2 
| ——— E/(10°GPa) -7 = w(x) (5) 
--- al(10°K") a? nog és 
4b $ we where u(x) and wọ (x) are the axial displacement and deflection on the 
a «/(10°W/mK) ae geometrical middle plane, respectively. p(x) is the rotation angle about 
Q the x-axes. (z) denotes the shape function determining the distribution 
5 of transverse shear strain and stress along the thickness. We have its 
a 3 form in SSDT [27], the first-order shear deformation theory (FSDT) and 
Z the classical shear deformation theory (CSDT) as 
© 
a 5 sin | 2] ssor 
o m h 
E eee SS = P=) z FSDT 
A 0 CSDT (6) 
1 a ENE AENEA a Sse ae Under the framework of small deformation, the strain components 
Be et ee can be expressed as 
i Ex = Uy’ (x) — Z — Cog Œ) + Y) — Clp Œ), Ve = PO E) 
300 (7) 
TK So the stress components in the k layer are obtained as 
i i i i i (k) 
Fig. 2. Material properties of constituents with respect to temperature change. o®) = BWe,,, 20 = x z te a = —ROQWMAT m 
+9 
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where Co is equal to 0 for homogeneous isotropic beams. As shown in 
Fig. 3, it is worth noting that Cp and C; are also temperature-dependent, 
so the physical neutral plane would vary with thermal environments. 
The displacements ux and u, of the beam in the x and z directions could 
be given by 
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Fig. 3. Flowchart of iterative solution for TD case. 
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where AT = T — h and ùh is the initial temperature. 
3.2. Buckling analysis 


The strain energy U and the potential energy V due to thermal stress 
can be expressed as 


3 
1 L hk 
U=>b SES _, Cv Eee + HE Ye) ded a 
k=1 


3 
1, pl he O p 
v=, 2AA x w, (x)dzdx 
k=1 ~ 


Applying Hamilton principle ô T= 6(U + V) = 0 yields 


(10) 
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For clamped-clamped sandwich FGM beams, the boundary condi- 
tions are 


Uo (X) = w(x) = w’ (x) = Y(xX)=0, atx=0,L 


Here, to satisfy the boundary conditions, the displacement compo- 
nents are assumed as [2,4]. 


U(x) =} &,sin (Ax)cos (Ax) 
w(x) =r sin (Ax)sin (Ax) 


Po) =n- sin (Ax)cos (Ax) (12) 
where 4 = mz/L, and &,(j = u, œ, p) are the unknown amplitudes. 

Substituting Eq. (12) into Eq. (11), the generalized eigenvalue 
problem for thermal buckling can be written as 


{[K] + [KT] }{A} = 0 (13) 


where {A} = {ép En» cae and [K] and [KT] can be defined as 
4RAn — 8B 40°C 0 0 0 

[K] = | — 84B 1644D, -8F, |, [K= fo 422NT | 
4VCy — 8RFy 4?Hi + Ass 0o 0 0 


h 
fAu, Bu, Du}= Da Sa t EV{1, z— Co, (z — Co)?}dz 


h 
{Ci Ar A= dy, E k E® 


WE) -= Cr E- COPE) - GI, WO - Paz 
Ass= Dee fe SE i E [v (z) Pdz 
NT= Dear Spe, -Ea (T — dz 
(14) 


Additionally, it is worth mentioning that the stretching-bending 
coupling effect would disappear after introducing the physical neutral 
plane concept, which leads to By, = Cu = 0 in Eq. (14), so that the 
calculation can be simplified. 


4. Thermal environments 


Assume the FGM beams are exposed to thermal environments, in- 
cluding uniform temperature rise (UTR), linear temperature rise (LTR) 
and nonlinear temperature rise (NLTR). 

For UTR, the temperature rises AT from the initial temperature 
To = 300K to the temperature T, so the temperature field can be ex- 
pressed as 


T= + AT (15) 


For LTR, the temperature field distributes linearly from the top 
surface T, to the bottom surface Tm = 300K, through the thickness, and 
it can be expressed as 

+5) 
h 


For NLTR, the temperature field distributes nonlinearly from the top 
surface T, to the bottom surface Tm = 300K. In this case, the one-di- 
mensional steady-state heat conduction through the thickness is con- 
sidered, i.e. 


_d mo) zoj- 
ak (z) 0 


T(2) =T™+(h- (+ 


(16) 


(17) 


Taking into account the thermal continuous conditions between 
different layers yields 
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In the case of FGM-I, Eq. (18) can be solved by using an approx- 
imation of polynomial series expansion. 
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where r denotes the item numbers of the series and r = 5 is chosen to 
satisfy the calculation precision. 

To examine the temperature dependence effect of material proper- 
ties on thermal buckling of porous sandwich beams, a comparison will 
be made between the temperature-dependence (TD) and temperature- 
independence (TID) cases. For TD case, the thermal-mechanical cou- 
pling buckling solution is solved by the iterative procedure as shown in 
Fig. 4. The determinant of the coefficient matrix in Eq. (13) must vanish 
to ensure non-zero solution, giving [K] + [KT] = 0. Specially, the cri- 
tical buckling temperature AT., = T — % under UTR can be solved 
analytically as 


_ 1644(F — Dy Hy) — M?Ass Dn 


(20) 


where 


3 
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Solve AT!" by Eq.(13) with 
the properties at T = T, + AT} 


i+] i 
ATs" -AT OoN 


cr 


Ala = AT," 


Fig. 4. Flowchart of iterative solution for TD case. 


5. Theoretical verification and result analyses 
5.1. Comparison study 


To verify the present theory, comparisons are made for perfect FGM 
beams and perfect sandwich FGM beams subjected to UTR thermal 
environment. Setting h, = 0 would degenerate the sandwich beam into 
a single FGM layer. From Table 2, the comparison study indicates the 
critical buckling temperatures, predicted by the present formulation 
based on physical neutral plane, are in agreement with the theoretical 
results of Ref. [4,28] under both TID and TD cases and the deviations 
are within 0.1%. For perfect sandwich FGM beams without porosity, an 
ABAQUS finite element model was built to verify the present theory. As 
the aforementioned, FGMs can be seen as laminated structures with 
multiple composite layers. For modeling such a laminated structure, the 
Composite Layups tools in ABAQUS is used to create a conventional 
shell with 50 layers. Each layer is homogeneous and has the same 
thickness, but the material properties of the laminated structure vary 
spatially in a gradient. Here the element type is S8R and the Lanczos 
method in the linear perturbation buckle step is used to solve the cri- 
tical temperatures. Table 3 indicates a maximum deviation of 0.3% 
between the present theoretical predictions and those by ABAQUS. 
Tables 2 and 3 also reveal that the critical temperatures predicted in the 


Table 2 
Verification of the critical temperature (in K) of FGM beams (L/h = 25, $ = 0, 
UTR). 


Cases Sources n=0.5 n=1 n=2 n=5 n=10 

TID Ref. [28] 509.890 458.680 423.530 394.390 376.140 
Ref. [4] 510.002 458.780 423.624 394.488 376.223 
Present (GMP) 510.189 458.904 423.629 394.393 376.191 
Present (PNP) 510.189 458.904 423.629 394.393 376.191 

TD Ref. [28] 399.500 367.320 345.150 326.400 313.580 
Ref. [4] 399.571 367.388 345.218 326.462 313.648 
Present (GMP) 399.694 367.470 345.213 326.380 313.617 
Present (PNP) 399.694 367.470 345.213 326.380 313.617 
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Table 3 
Verification of the critical temperature (in K) of sandwich Siz3N4/SUS304 FGM 
beams (L/h = 50, TID). 


ha/he N ABAQUS Present 
SSDT, PNP SSDT, GMP 

0.125 0.2 137.070 137.121 137.121 

1 116.460 116.556 116.556 

5 106.120 106.205 106.205 
0.5 0.2 126.590 126.932 126.932 

1 117.350 117.447 117.447 

5 112.580 112.638 112.638 

Table 4 


Effects of slenderness ratio on the critical temperature (in K) of sandwich FGM-I 
(n = 1) beams under different thermal environments (TID). 


Lih B ha/hş = 0.125 ha/hş = 0.5 
SSDT FSDT CSDT SSDT FSDT CSDT 

UTR 

30 0 321.520 321.577 325.045 323.997 324.049 327.516 
0.1 360.756 360.860 364.861 354.512 354.685 358.622 
0.2 407.968 408.149 412.819 387.439 387.792 392.268 

50 o0 116.556 116.564 117.016 117.447 117.453 117.906 
0.1 130.814 130.828 131.350 128.567 128.590 129.104 
0.2 147.981 148.005 148.615 140.586 140.632 141.217 

LTR 

30 0 678.490 678.610 685.930 690.703 690.814 698.206 
0.1 766.466 766.687 775.188 760.309 760.679 769.124 
0.2 873.955 874.343 873.955 836.358 837.119 836.358 

50 0 245.964 245.980 246.935 250.375 250.390 251.354 
0.1 277.929 277.958 279.068 275.734 275.782 276.885 
0.2 317.008 317.059 317.008 303.480 303.580 303.480 

NLTR 

30 0 712.988 713.114 720.806 732.405 732.523 740.361 
0.1 808.863 809.097 818.068 809.345 809.739 818.728 
0.2 926.767 927.178 937.787 894.299 895.113 905.446 

50 0 258.470 258.486 259.490 265.492 265.507 266.530 
0.1 293.303 293.334 294.504 293.517 293.569 294.742 
0.2 336.164 336.219 337.603 324.504 324.612 325.961 


geometrical middle plane (GMP) and the physical neutral plane (PNP) 
cases are nearly equal to each other, which is consistent with the results 
in Ref. [17]. Therefore, the following theoretical results will be pro- 
duced by the formulation based on PNP. 

From Table 4, the critical temperatures of Si3N4/SUS304 FGM-I 
sandwich beams predicted by SSDT and FSDT (the shear correct factor 
is mata) are slightly smaller than those of CSDT. The deviations are 
within 1.2%. Meanwhile, with the increasing value of the h,/h; and 
L/h, the critical temperature Ter decreases greatly. However, it is in- 
teresting that T, would rise as the porosity parameter increases. For 
instances, when the porosity increases by 0.1, the critical temperature 
increases by 9-15% for the beams with h,/h~ = 0.125. 


5.2. Thermal buckling under UTR 


In the following parts, SSDT is used to calculate the theoretical re- 
sults and all the calculation parameters are given in the figures for 
brevity. Fig. 5 gives the changing rules of the critical temperature with 
the power-law index and porosity parameter. Fig. 5 (a) and (b) are 
plotted for SisN4/SUS304 and ZrO2/SUS304 FGMs, respectively. T., 
decreases monotonically with the rising value of n for Si3N,/SUS304 
FGMs, but increases for ZrO2/SUS304 FGMs. Meanwhile, these curves 
would tend to be a steady value in a large value of n. It means that, for 
ZrO2/SUS304 FGMs, the existence of the ceramic constituent would 
reduce the critical temperature, rather than increasing the thermal re- 
sistance of structures. So it is confirmed that to improve the critical 


Y. Liu et al. Composites Part B 168 (2019) 236-242 
(a) 220 (b) 120 
200) & Si,N,/SUS304 110 ZrO,/SUS304 L/h=50, FGM-I, 
UTR, h,/h=0.125 
180 
-->- B=0,TID 
aie —— f=0,TD 
K 
S 140 --+- £=0.1,TID 
120 —— f=0.1,TD 
iod -->- B=0.2,TID 
—— f=0.2,TD 
80 1 1 4 4 4 1 50 1 4 1 4 1 1 
0 1 2 3 4 5 0 1 2 3 4 5 
n n 
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Fig. 6. Effects of the porosity and power-law index on thermal buckling of sandwich FGM-I beams under LTR and NLTR. 


temperature, the ceramic constituent ZrO% is unsuitable to mix with 
SUS304. This is caused by its large thermal expansion coefficient re- 
lative to other ceramic constituents as witnessed in Fig. 2. Therefore, to 
improve the thermal buckling load of FGM beams, ceramic constituents 
with lower thermal expansion coefficients would be preferred. 

Besides, the TID results are approximately 6-25% higher than those 
of TD, which indicates that neglecting the temperature-dependence of 
the constituent would overestimate the actual critical temperature. 
Meanwhile, from both figures, it is interesting that the existence of the 
porosity would heighten the overall curves of critical buckling tem- 
peratures. This result is quite different from the changing rule of the 
critical buckling mechanical load of porous FGM plates obtained in 
Refs. [12,13]. In fact, it is not paradoxical because porosity would de- 
crease the Young's modulus and the thermal expansion coefficient of 
FGMs synchronously. Reduction in the latter would certainly decrease 
the load induced by temperature rise before thermal buckling. 


241 


5.3. Thermal buckling analyses under LTR and NLTR 


Practically, FGMs are usually served for thermal barrier environ- 
ments with temperature gradients, where temperature-dependent ma- 
terial properties of FGMs should be considered as the aforementioned. 
To focus on the temperature gradient effects, both LTR and NLTR are 
discussed in this part. Fig. 6 illustrates the critical temperature chan- 
ging with the porosity parameter and the power-law index of Si3N4/ 
SUS304 FGMs with TD material properties. Similar to Fig. 5(a), To 
decreases monotonically with the rising value of n. 

It is worth noting that there are significant differences between the 
LTR and the NLTR case, and the latter includes the practical thermal 
conduction effect. The NLTR results are higher by 4-17% than those of 
LTR. This illustrates that, as a simplified method of thermal environ- 
ment, LTR would arouse tremendous error and greatly underestimate 
the critical buckling temperature. Besides, comparing Fig. 6 (a) and (b) 
indicates that the critical temperature of the sandwich beam with a 
smaller face-to-core ratio h,/h¢ is more sensitive to the change of the 
power-law index and the porosity parameter. 
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Fig. 7. Porosity effects on thermal buckling of sandwich FGM beams under 
NLTR. 


5.4. Porosity effects 


To further reveal the effects of porosity on thermal buckling of 
sandwich FGM beams in the practical conditions, by considering the TD 
and NLTR cases, comparisons are made between different porosity 
distributions including the FGM-I and FGM-II cases as shown in Fig. 7. 
In both cases, T., increases with the rising value of the porosity para- 
meter, which is similar to the aforementioned. The critical buckling 
temperatures in the FGM-I case are larger than those in the FGM-II case 
and the deviations enlarge gradually as the porosity increases. Hence, 
in the FGM-I case, porosity effects more significantly on the critical 
temperatures of FGM sandwich beams than in the FGM-II case. 


6. Conclusion 


Thermal-mechanical coupling buckling behaviors of the clamped- 
clamped porous FGM sandwich beam is investigated in this paper, from 
which several conclusions can be drawn. 


1. The critical temperature of FGM sandwich beams predicted by the 
theory based on the physical neutral plane is nearly equal to that 
based on the geometrical middle plane. 

. Neglecting the temperature-dependence of the constituent materials 
would overestimate the actual critical temperature. Meanwhile, an 
approximation of the actual thermal environment by a linear tem- 
perature rise method would arouse tremendous error and greatly 
underestimate the prediction of the critical buckling load. Hence, 
thermal buckling analysis for FGM sandwich beams should take the 
temperature-dependent material properties and the practical 
thermal conduction into consideration. 

. The existence of porosity would increase the critical temperatures of 
FGM sandwich beams, and the uniform porosity distribution (FGM- 
I) affects more significantly on the increase of critical temperatures 
than the non-uniform porosity distribution (FGM-ID. 

. The critical temperature of the sandwich beam with a smaller face- 
to-core ratio is more sensitive to the change in the power-law index 
and the porosity parameter. 

. To improve the thermal buckling load of FGM beams, the ceramic 
constituent with a lower thermal expansion coefficient would be 
preferred. 
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